Spin Wave Excitations in AFei. 5 Se 2 (A=K, Tl): Analytical Study 
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We have analytically solved the spin wave excitations for the intercalated ternary iron-selenide 
AFei.5Se2 (A=K, Tl) in the 4x2 collinear antiferromagnetic order. It is found that there are one 
acoustic branch (gapless Goldstone mode) and two gapful optical branches of spin wave excitations 
with each in double degeneracy. By examining the non-imaginary excitation frequency condition, 
we can determine the corresponding phase boundary. The exchange couplings between Fe moments 
in AFei.5Se2 are derived based on the first-principles total energy calculations. The Fe spin is found 
to be S — | through computing the antiferromagnetic quantum fluctuation. And it is further 
found that a very small spin-orientation anisotropy can remarkably suppress the antiferromagnetic 
quantum fluctuation. The spin dynamical structure factors are calculated and discussed in associated 
with neutron inelastic scattering experiment. 



PACS numbers: 74.70.Xa, 75.30.Ds, 75.30.Et 



I. INTRODUCTION 

The discovery of iron-based superconductors 1 has 
stimulated great interest on the investigation of uncon- 
ventional superconducting mechanism, in which mag- 
netism is considered to play a substantial role. It is 
well-known that the parent compound of a cuprate super- 
conductor is an antiferromagnetic (AFM) Mott insulator 
with a checkerboard AFM order on the copper square lat- 
tice. In contrast, the parent compound of an iron-based 
superconductor was found to be an AFM semi-metal 2 
with either collinear— ^ or bi-collinear— ~— AFM order on 
the iron square lattice. Regarding the nature of the mag- 
netism, there are basically two contradictive views. The 
one& is based on itinerant electron picture, in which the 
Fermi surface nesting is responsible for the collinear AFM 
order. On the contrary, the other one is based on local 
moment interactions which can be described by the J\- 
J<i frustrated Heisenberg modeli^^ And it was further 
showni that the underlying driving force herein is the 
As-bridged AFM superexchange interaction between a 
pair of next-nearest-neighboring fluctuating Fe local mo- 
ments embedded in itinerant electrons. There are now 
more and more evidences in favor of the fluctuating Fe 
local moment picture. Especially, the neutron inelastic 
scattering experiments have shown that the low-energy 
magnetic excitations can be well described by the spin 
waves based on the J\-Ji Heisenberg modelAi~— 

The newly discovered intercalated ternary iron- 
selenide superconductors A 2/ Fe a; Se2 (A=K, Tl)^£r— show 
rich phase diagrams and many unusual physical proper- 
ties that have not been found in other iron-based su- 
perconductors, for example, the superconductivity was 
found to coexist with a strong AFM order with a giant 
magnetic moment of 3.31 /is/Fe formed below a Neel 
temperature of 559i f 18 i 19 and to be proximity to an AFM 
insulating phased These reveal the close relationship 
between unconventional superconductivity and antiferro- 
magnetism once more, and have triggered another surge 



of interest for the investigation of unconventional super- 
conducting mechanism. 

The compounds A a Fe x Se2 have the ThCr2Si2 type 
crystal structure, isostructural with 122-type iron 
pnictidesi2£ However, the stable structures of A y Fe ;E Se2 
contain Fe vacancies ordered in either \fh x v5 or 4 x 2 
superstructure due to the balance required in chemical 
valences, which respectively correspond to A0.sFe1.6Se2 
with one- fifth Fe vacancies or AFei.5Se2 with one-quarter 
Fe vacancies. For A0.sFe1.6Se2, the neutron observation 
has found that it has a \J% x \fh blocked checkerboard 
AFM order to match the Fe vacancy superstructure^ 
For AFei.5Se2, the first-principles electronic structure 
calculations predicted 21 that its ground state is in a 4 x 2 
collinear AFM order, as shown in Fig. Hid). More- 
over, the calculations further showe d 21 ' 22 that both com- 
pounds AFei.5Se2 and A0.sFe1.6Se2 are antiferromagnetic 
semiconductors with band gaps of dozens and hundreds 
meV, respectively. Such band gaps have been confirmed 
by the recent ARPES and transport measurements^ It 
was further proposed that the parent compound of an 
Aj / Fe a; Se2 superconductor is an AFM semiconductor ei- 
ther A0.8Fe1.6Se2 or AFei.5Se2i 21 ' 22 The latest neutron 
diffraction experiment shows that the parent compound 
is likely AFei. 5 Se 2 with a 4 x 2 collinear AFM order ,21 In 
order to well understand the magnetism in the compound 
AFei.5Se2, we have studied the magnetic excitations and 
spin dynamical structure factors, which can be directly 
detected by neutron experiments. Our approach is based 
on the linearized spin wave theory upon the effective spin 
Heisenberg model. 



II. EFFECTIVE MODEL 

As schematically shown in Fig. [IJd), in the ground 
state of the compound AFei.sSe2r^ the Fe vacancies are 
ordered in a rhombus structure, in which unit cell there 
are two inequivalent Fe atoms according to the number 
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of neighboring Fe atoms, namely 2-Fe-neighbored and 3- 
Fe-neighbored Fe atoms respectively. The corresponding 
magnetic order is a 4 x 2 collinear AFM order (also called 
A-col AFM order here), as shown in Fig. [TJd), in which 
the Fe moments are antifcrromagnetically ordered along 
the lines without Fe vacancies and ferromagnetically or- 
dered along the lines perpendicular. 

To describe the magnetism in the compound 
AFei.5Se2, considering the Se-bridging effect^! we 
adopt the spin Hcisenberg model with the nearest and 
next-nearest neighbor exchange couplings between the Fe 
moments with quantum spin S as follows, 

H = Ji Si ■ Ss a + Jib J]] Si ■ Ss b + J2 Si ■ Sj, 

<i,<5a) (i,S b ) <ii» 

whereas (i,5 a ), (i,5b), and <C ij 3> denote the summa- 
tion over the nearest neighbors along the lines without Fe 
vacancies, the nearest neighbors along the lines with Fe 
vacancies, and the next-nearest neighbors, respectively. 
To be more specific, the exchange coupling J\ a links a 
pair of 3-Fe-neighbored Fe atoms and the exchange cou- 
pling Jib links a 2-Fe-neighbored Fe atom and a 3-Fe- 
neighbored Fe atom, while the exchange coupling J 2 con- 
nects a 2-Fe-neighbored Fe atom with a 3-Fe-neighbored 
Fe atom (see Fig. Hfd)). 

The magnetic phase diagram of the Hamiltonian (JTJ 
has been studied by classical Monte Carlo simulations 
in Ref. [3ll In the shadowed part as shown in Fig. [2j 
the ground state of the Hamiltonian (JTJ is in the A-col 
AFM phase with a 4 x 2 collinear AFM order, whose spin 
dynamics will be studied below. 



III. SPIN DYNAMICS 

For a spin Heisenberg model with a long-range 
magnetic order in its ground state, the linearized 
spin wave theory with the Holstein-Primakoff (HP) 
transformation 25 is a standard approach to obtain the 
spin wave excitations and other dynamical properties. In 
the case of a simple magnetic unit cell, namely no more 
than two spins per cell, the approach can easily and di- 
rectly give an analytical solution, for example, the well- 
known quadratic and linear dispersion behaviors in low 
energy for ferromagnetic and antiferromagnetic orders, 
respectively. However, for a complex magnetic struc- 
ture with a magnetic unit cell containing more than two 
spins, it is still a severely challenging task to analytically 
solve the spin wave excitations even though lots of ef- 
forts have been devoted in the past One thus has 
to be satisfied with numerical solutions. The underlying 
difficulty is as follows. After the linearized HP transfor- 
mation, the spin Heisenberg model is transformed into a 
quadratic Bosonic Hamiltonian. Conventionally one at- 
tempts to construct a Bogoliubov transformation to di- 
agonalize the Bosonic Hamiltonian. 3 - In general, there 
is no practical procedure to analytically construct such 
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FIG. 1. (Color online) Different magnetic orders on an Fe- 
Fe square layer with one-quarter Fe vacancies ordered in 
rhombus: (a) ferromagnetic order; (b) Neel antiferromagnetic 
(AFM) order, in which all pairs of the nearest-neighboring Fe 
moments are in antiparallel order; (c) P-collinear AFM or- 
der, in which the Fe moments are in antiparallel order along 
the lines with the Fe vacancies; (d) A-collinear AFM order 
(ground state), in which the Fe moments are in antiparallel 
order along the lines without the Fe vacancies. The filled 
circles represent the Fe atoms while the empty circles rep- 
resent the Fe vacancies. The red (blue) arrows represent the 
up-spins (down-spins). The A-collinear AFM order can be di- 
vided into six ferromagnetic sublattices, labeled successively 
as ai to a?,. J\ a and Ju are the exchange couplings along 
the nearest neighbor Fe-Fe directions without and with the 
Fe vacancies, respectively. J2 is the next-nearest neighbor ex- 
change coupling. The rectangle enclosed by the thick blue 
solid lines denote a 4 x 2 magnetic unit cell. The x and y axes 
are also shown. 



a Bogoliubov transformation for a Bosonic Hamiltonian 
with more than two component bosons, corresponding to 
a case of more than two spins per unit cell. 

Nevertheless an analytical solution is usually in desire 
and very helpful to understand physics, especially more 
helpful in the case of a spin Heisenberg model in com- 
parison with neutron measurement, for example, to de- 
termine the exchange couplings. Here by studying the 
Hamiltonian (fTJ), rather than constructing a Bogoliubov 
transformation as conventionally, we are going to illus- 
trate an alternative method, which is able to analytically 
solve spin wave excitations for the cases of multi-spin 
unit cells. This method employs the equation of motion 
to construct a secular equation, i.e. an algebraic equa- 
tion. According to the algebraic basic theorem, an alge- 
braic equation can be analytically solved up to the fourth 
power, which corresponds to an antiferromagnetic order 
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with eight spins per unit cell. To our knowledge, the ad- 
vantage of this method applied on a complex magnetic 
structure has not been well recognized. 
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FIG. 2. (Color online) Phase diagram of A-collinear antiferro- 
magnetic order, marked by the shadowed region, in variance 
with the exchange couplings Ji a , Jit, and J2, as defined in 
Fig. 0Jd) • Black solid and red dashed lines represent the clas- 
sical Monte Carlo calculation (Ref. I3H ) and our result deter- 
mined by examining the non-imaginary excitation frequency 
condition, respectively. Blue and green solid dots stand for 
the positions of KFei.sSe2 and TlFei.sSe2 respectively. 



A. Spin wave excitations 

We now consider the Hamiltonian ([T]) in the A-col 
AFM phase with a 4 x 2 collinear AFM order. As shown 
in Fig. [ljd), the magnetic unit cell contains two vacan- 
cies and six spins, i.e. three up spins and three down 
spins, which are labeled by with £ (£ = 1, . . . , 6) be- 
ing odd number or even number, corresponding to up or 
down spins respectively. In other words, the original an- 
tiferromagnetic lattice can be divided into two vacancy- 
sublattices and six ferromagnetic sublattices, in each of 
which we can perform the linearized Holstein-Primakoff 
transformation for the quantum spin S at each site i as 
follows, 




S? = S- 

£ G odd number 



where (ctji) is the £-th component boson creation (an- 
nihilation) operator and i is the site index in the £-th 
sublattice. We then perform the following Fourier trans- 
formation, 
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where denotes the position vector of site i in the £-th 
sublattice, k = (k x , k y ) is a wave vector in the magnetic 
Brillouin zone, and N is the total site number of each 
sublattice. 

It turns out that the spin Hamiltonian ([1} is now trans- 
formed into a quadratic Bosonic Hamiltonian in the mo- 
mentum space as follows^ 



H = 

E 

#1 



Eo + H 1 , 
= (4Ji 6 - 4Ji Q 



8J 2 )NS(S + l) 



(2) 



whereas 
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where 2 is the length of two nearest neighbor sites, namely 
the Fe-Fe bond length. In the following we take I as the 
length unit for convenience. 

Next we construct a secular equation to solve the eigen- 
values of the Hamiltonian ([2]) by using the equation of 
motion method. We first consider the following linear 
combination, 



£Godd 



i k? a| k 



CGodd 



^2 u k?a£k, 

even 



w ^4k- 

££even 



Without loss of generality, the coefficients w k j and vu^ 
are assumed to be complex numbers. Here d k and /3 k 
are the boson operators, which satisfy the relationship 



of the boson commutation, namely [dk,dL] 



?kk' 



and 



[/3k> /? k /] = <^ k k'- This imposes the following constraints 
upon the coefficients u k j and u k( r, 

- J2ieodd U H U U + Eleven U H U U = X ' 
Ef^odd "k?^ - Eleven = 1, 

whereas * means the complex conjugate and remains the 
same meaning as follows. 

Now we assume that d k or /3 k describes an eigen-mode 
of the Hamiltonian ([5]), it is then required that d k or /3 k 
fulfils the following equations of motion respectively, 

ifr&k = [dk, H) = A k a k or £7i/3 k = [/3k, #] = A k /3 k , 
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where Ak and A k are the respective eigenvalues. This 
gives rise to a generalized eigenvalue problem at each 
wave vector k in the magnetic Brillouin zone, 

H^U k = -A k cr 3 f7k, U k = (u k i,u k 2,u k3 ,u k4: ,u k5 ,u k6 ) T , 

(3) 

or 

H%*Vk = A k Cr 3 Vk, Vk = (Wkl,W k 2,Wk3,'fk4,Wk5,'yk6) T , 

(4) 

where the 6x6 matrix is given in Eq. ([2|) , and 03 is 
a 6 x 6 diagonal matrix with diagonal elements being (1, 
-1, 1, -1, 1, -1) respectively. 

Eqs. (0) and dU will show that each eigenvalue of the 
Hamiltonian ([2]) is in double degeneracy, which reflects 
the symmetry or equivalence between the spin-up and 
spin-down in the ground state. For convenience, we only 
deal with Eq. ([3]) in detail. We will not diagonalize Eq. 
d3j directly since this yields a set of quadratic equations 
with at least 36 unknowns, which is unlikely to give an 
analytical expression. Instead, we come to solve the sec- 
ular equation of the Hamiltonian ([3]), namely zeroing the 
determinant of the matrix + Ak<7 3 , which is the con- 
dition for there being a nonzero vector U k . Actually, the 
Hamiltonian (jj) has the same secular equation. The sec- 
ular equation turns out to be an algebraic equation with 
one unknown Ak, namely 

A k + b\t + c\l + d = 0, (5) 

whereas the coefficients 6, c, and d all are real functions 
of the wave vector k, and composed of the exchange cou- 
plings. For the detailed expressions on 6, c, and d, please 
refer to Appendix I. 

Setting 7 = A k , we then have 7 3 + bj 2 + cj + d = 
from Eq. ([5]). Further setting 7 = a — |, the equation 
reduces to a 3 +pa + q — 0, whereas p = — y + c, and q = 

ff-y+d. When the discriminant A = [(|) 2 +(|) 3 ] < 0, 
there are three real roots for this cubic equation, namely 
a% = 2-^rcos(p, a?. = 2tyr cos(<fi+=§-), 013 — 2yfr cos(<^>+ 
whereas r — (|) 3 and <j> = ^ arccos(— ^). Now 
we arrive at the condition, which makes all six roots of 
Eq. §5§ being real, is that a\ — | > and a-i — | > 
and «3 — I > 0. In the end, we analytically obtain six 
branches of eigenvalues for the Hamiltonian ([2]) , in three 
pairs of positive and negative, similar to the case of calcu- 
lating lattice phonon excitations. And we keep the three 
positive branches to describe the spin wave excitations 
as follows, 

Aik = ?kt>i(k) = Sy/a.\ — 6/3, 
A 2k = foJa(k) = Sy /al - 6/3 , (6) 
A 3k = ficj 3 (k) = S yjaz - 6/3. 

Likewise we can analytically solve Eq. (J4j) to obtain 
the same set of eigenvalues, but with the different eigen- 
vectors. Lets take the complex conjugate of Eq. ((3]), 
namely H**U£ = -A k r7 3 L/ k . Comparing with Eq. ([J}, 



we can see that a positive or negative eigenvalue (±A k ) 
with eigenvector £/ k of Eq. © is a negative or positive 
eigenvalue (=FA k ) with eigenvector Vk = U£ of Eq. (jH). 
This directly shows that the spin wave excitations are in 
double degeneracy. 

The Hamiltonian © is derived from the spin Hamil- 
tonian (JT|) on the assumption of the ground state being 
in the 4 x 2 collinear AFM order (see Fig. [Tfd)). As 
shown above, the eigenvalues of the Hamiltonian ([2} may 
be imaginary. Such a case, if happening, indicates that 
the 4x2 collinear AFM order is unstable for the spin 
Hamiltonian (JTJ). This is similar to the case of imagi- 
nary phonon vibration modes, which indicates that the 
corresponding crystal structure is unstable. Thus we can 
determine the boundary of A-col phase in the phase dia- 
gram of the spin Hamiltonian (JTJ) by checking the condi- 
tion whether the eigenvalues of the Hamiltonian @ being 
real or not. Figure [5] shows our calculated phase bound- 
ary by the red dashed line, which is quite close to the 
one given by the classical Monte Carlo simulations. The 
difference between the two is that the AFM quantum 
fluctuations are included in our spin wave calculations 
but not in the classical Monte Carlo simulations. 

In order to calculate physical quantities, we need to 
know the eigenstates corresponding to the spin wave ex- 
citations, which but are not directly U k and V k of Eqs. 
([3]) and ([J} because of the boson's characteristic. Let 
us denote the eigenvectors corresponding to the positive 
and negative eigenvalues of Eq. ([3]) by and K k 

respectively, with 77 = 1,2,3 being the branch index of 
the spin wave excitations. Then the Hamiltonian ^ will 
be diagonalized into 

H = E Q + J2 ^(ajk^k + )8{iAk + 1) 

v/k 

after the following linear transformation, 

<V = - E 5eodd «S<4k + Eleven «$Ak> 
&rik = E^eodd M ij£,k^k — Eleven U ij£,k^k- 

The above set of equations can be rewritten in a matrix 
form as 

(«ik /3 lk "2k /?|k "3k /3;j k ) = (a\ k a 2k 63k a4k a\ k a 6 k)Qk- 

(7) 

Here the six columns of the matrix Q k are nothing but 
the eigenstates corresponding to the spin wave excita- 
tions. 



B. Spin dynamical structure factor 

As shown above, S z just contributes to elastic scat- 
tering rather than inelastic scattering since it does not 
change the number of magnons. In contrast, S x and S y 
contribute to inelastic scattering by changing the num- 
ber of magnons. In zero temperature the spin dynamical 
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structure factor (SDSF) in inelastic scattering process 
through single magnon excitations is defined as 

5(k,«) = E/E—K/I^C^lo)! 2 ^-^) 

= sE 7 IQ~kl 2 ^-^/), [ j 

where |0) is the vacuum state, |/) denotes the final states 
of a spin system with excitation energy Wf, and Q~\ 
means the sum of the elements in the 7-th row of the 
matrix Q^ 1 given in Eq. ([7]). 



TABLE I. Calculated energies of the four different magnetic 
states, i.e. the ferromagnetic state (Ef), Neel antiferromag- 
netic state (En), P-collinear antiferromagnetic state (Ep), 
and A-collinear antiferromagnetic state (Ea), as shown in Fig. 
[TJ for AFei. 5 Se 2 (A=K, Tl) (unit: meV/Fe). The energy of 
the nonmagnetic state is set to zero. The calculated exchange 
couplings Ji a , Jii,, and J2 are also listed (unit: meV/S 2 ). 



Compounds 


E F 


En 


E P 


E A 


Jia 


Jib 


J 2 


KFei. 5 Se 2 


-157.1 


-253.0 


-325.3 


-370.0 


52.7 


19.2 


53.5 


TIFei.sSea 


-79.4 


-176.0 


-255.7 


-283.8 


46.8 


25.7 


53.3 



C. Sublattice magnetization 

In the £-th sublattice the AFM quantum fluctuation re- 
duces the staggered magnetization from its classical value 
S by the following quantity AS^, 

A <% = S- ^(EtilS-a^M) = ^(Ek4k a «k) 
= |fEk E w =i,3,5 IQ^kl 2 ' for i e odd number, 
I Tf Ek E^=2,4, 6 IQ^.kl 2 ' for £ e even number. 

(9) 

where denotes the element of the matrix Q k 1 at 

the crossing of the i/-th row and the £-th column. The 
summation of k can be replaced by a two dimensional 
integral over the whole Brillouin Zone. It is well-known 
that the AFM quantum fluctuation is strongly influ- 
enced by the spatial dimension and the number of near- 
est neighbors of a site, namely the coordination number. 
There are two kinds of two Fe atoms (sites) in a mag- 
netic unit cell in AFei.5Se2, one is 2-Fe-neighbored and 
the other is 3-Fe-neighbored. Thus we need to calcu- 
late the AFM quantum fluctuation on these two kinds 
of sites separately. Moreover, there usually exists small 
spin-orientation anisotropy in realistic materials. To ac- 
count for such anisotropy, one standard approach is to 
add the term of —A z E Sf 2 into the Hamiltonian ((T|). 
This turns out to add 2A 2 to all diagonal term of H\ 
and -3A Z NS(S + 1) to E in Eq. ^ respectively. 



IV. RESULTS AND DISCUSSION 

The above analytical study will be very helpful both 
experimentally and theoretically. Experimentally in com- 
parison with neutron inelastic scattering on AFei.5Se2, 
we can determine the exchange couplings Jy a , J\ , and 
J2 through the spin wave excitations (Eqs. ((6]) and (0). 
This will help us better understand the magnetism in 
AFei.5Se2- On the other hand, these exchange couplings 
can be theoretically derived from the relative energies 
of different magnetic states with respect to the non- 
magnetic state. For the detailed derivation, please refer 
to the appendix in Ref. 0- To be specific, we need to first 
obtain total energies of the four different magnetic states, 
i.e. the ferromagnetic state (Ep), Neel AFM state (En), 
P-collinear AFM state (E P ), and A-collinear AFM state 



(Ea), as shown respectively in Fig. [TJ The energy differ- 
ences among these different states result in three linearly 
independent equations as follows, 

E N -E F = -4( Jia + Jib)/3, 
E A -E N = 4( J lb - 2 J 2 )/3, 
Ea Ep — -4(J lo - J lb )/3, 

from which the exchange couplings J\ a , Jib, and J 2 can 
be uniquely determined. 

To obtain the total energies, we have carried out the 
first-principles electronic structure calculations on the 
compounds AFei.sSe2, which were reported in Ref. [2l|. 
The calculated total energies and subsequently derived 
exchange couplings for the compounds AFei.sSe2 are 
listed in Table fl] Accordingly, in Fig. [5] we mark the 
positions of AFei.sSe 2 in the magnetic phase diagram, 
both of which are in the A-col phase, but close to the 
phase boundary. Thus the ground states of AFei.sSe2 
both are in the 4x2 collinear AFM order, as observed in 
the neutron elastic scattering experiment, but they may 
be readily destructed by doping or applying pressure. 

From Table |H we also see that the exchange cou- 
plings J 2 are almost the same for both KFei.5Se2 and 
TlFei.5Se2, but the former is more anisotropic than the 
latter on the exchange couplings Jy, namely the differ- 
ence between Ji a and Ji . Such a magnetic anisotropy 
will help lessen the AFM quantum fluctuations. With the 
help of Eq. © , we are now able to calculate the AFM 
quantum fluctuations in AFei.5Se 2 . We summarize the 
calculated results in Table [TTJ As we find, the amount of 
the AFM quantum fluctuation in KFei.sSe2 is less than 
the one in TlFei. 5 Se 2 by about 20%. From Table II, we 
further find that the AFM quantum fluctuation on a 2- 
Fe-neighbored Fe atom is, as expected, rather larger than 
the one on a 3-Fe-neighbored atom, by about 10% and 
20% for KFei. 5 Se 2 and TlFei. 5 Se 2 respectively. Moreover 
a very small spin-orientation anisotropy can remarkably 
reduce the AFM quantum fluctuations, as shown in the 
case of A 2 = 0.02 J2, in which the AFM quantum fluctu- 
ation is already reduced from 0.2521 (0.3162) to 0.1943 
(0.2337) for KFei. 5 Se 2 (TlFei. 5 Se 2 ) by approximately 
20%. Experimentally the neutron measurements have 
shown that there is a small spin-orientation anisotropy 
in KFei.5Se2 and the other iron pnictides^^ ' 19 ' 24 In ad- 
dition, the effective ordering moment in KFei.5Se2 was 
found to be 2.8/is per Fe atom^ Here setting S = |, 
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TABLE II. Antiferromagnetic quantum fluctuations for 
AFei.5Se2 (A=K, Tl) with and without spin anisotropy 
respectively. ASW and AS^n denote the antiferromag- 
netic quantum fluctuation on a 2-Fe-neighbored and 3-Fe- 
neighbored Fe atom respectively. There are two 2-Fe- 
neighbored and four 3-Fe-neighbored Fe atoms in a 4 x 2 
magnetic unit cell (see Fig. Hd)). AS denotes the average 
antiferromagnetic quantum fluctuation on an Fe atom. 



Compounds 


A z = 0.0 


A z = 0.02J 2 


AS 2n 


AS 3n 


AS 


AA' 2 „ 


A5 3 „ 


AS 


KFei. s Se 2 


0.2685 


0.2439 


0.2521 


0.2093 


0.1868 


0.1943 


TlFei. 5 Se 2 


0.3665 


0.2911 


0.3162 


0.2780 


0.2115 


0.2337 



our calculations show that the ordering moment is 2.61/is 
with the spin anisotropy of A z = 0.02 J2 for KFei.5Se2, in 
good agreement with the neutron observation. We thus 
arrive with S = § for each Fe atom in KFe! 5 Se 2 . In the 
following discussion, we set S = § for each lattice site. 




FIG. 3. (Color online) Three branches of the spin wave ex- 
citations in A-collinear antiferromagnetic phase with a 4 x 2 
magnetic unit cell (see Fig. Old)), (a) KFei.5Se 2 , in which 
Jia = 23.4, Ji b = 8.5, and J 2 = 23.8 meV (S = §). (b) 
TlFei. 5 Se 2 , in which J lo = 20.8, J lb = 11.4, and J 2 = 23.7 
meV (S = §). Note k x is of spin-antiparallel direction with- 
out the Fe vacancies aligned and k y is of spin-parallel direction 
with the Fe vacancies aligned. Here the Fe-Fe bond length I 
is taken as the length unit for convenience. 



For an AFM structure with a magnetic unit cell con- 
taining six spins, physically there are six branches of spin 
wave excitations with every two in degeneracy, namely 
three 2-degenerated branches, among which one is the 
gapless Goldstone mode and the other two are gapful op- 
tical modes. With Table fl] and Eqs. ©, we can plot 
the spin wave spectra in the extended Brillouin zone for 
the compound AFei. 5 Se 2 in Fig. [31 As we see, there 
are exactly one gapless Goldstone mode and two gapful 
optical modes, and each branch of spin wave excitation 



is much more dispersive along the spin-antiparallel di- 
rection than along the spin-parallel direction. This is 
attributed to the Fe vacancies hindering the spin wave 
propagations, which are located along the spin-parallel 
direction rather than the spin-antiparallel direction (see 
Fig. [ljd)). Furthermore, such hindering becomes more 
severe, also the spin wave diffraction becomes weaker, 
for the higher spin wave excitations because of the cor- 
responding wavelengthes being shorter. Especially, the 
highest optical mode is almost dispersionless along the 
spin-parallel direction. 




k x k x k x 



FIG. 4. (Color online) Constant-energy cuts of the spin 
dynamical structure factor (SDSF) for KFei.sSe 2 , in which 
Jia = 23.4, Jib = 8.5, and J 2 = 23.8 meV (S = §). (a)- 
(f) are the SDSF at the excitation energies of 20, 40, 60, 80, 
100, and 120 meV, respectively. k x and k y both vary from 
-7r to 7r. Note k x is of spin-antiparallel direction and k y is of 
spin-parallel direction. 



To compute the spin dynamical structure factor by 
using Eq. (jSJ, the function 6(u> — u)f) is smoothed by 
a Gaussian distribution function, which is centered at 
the excitation energy ujf with the half- width at the half- 
maximum (HWHM) assigned to O.2J2. We plot the spin 
dynamical structure factor at a constant cut energy, rang- 
ing from 20 to 120 meV for both compounds, as shown 
in Figs. 4 and 5 respectively. 

For KFei.5Se2, referring to Fig. 3, one can see that 
the diffraction peaks come from the acoustic (Goldstone 
mode) spin waves for the excitation energy u)f < 50 meV, 
the acoustic and the first optical spin waves for 50 meV 
< ujf < 100 meV, and the first and second optical spin 
waves for u>f > 100 meV, respectively. When ujf = 20 
meV, Figure S] shows that the strongest diffraction peaks 
are located at (±7r, 0). In addition there are two small 
half elliptical diffraction peaks near T point and four al- 
most invisible diffraction peaks at the four corners of the 
extended Brillouin zone. When ujf increases from 20 to 
80 meV, these ellipses extend larger. At 100 meV, the 
diffraction peaks become line-shaped. When ujf = 120 
meV, the diffraction peaks become helical lines along the 
k y directions, namely spin-parallel direction. These fea- 
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(d) 



i;;isaoo 



FIG. 5. (Color online) Constant-energy cuts of the spin 
dynamical structure factor (SDSF) for TlFei.sSe2, in which 
Jia = 20.8, Jib = 11.4, and J 2 = 23.7 meV (5 = §). (a)- 
(f) are the SDSF at the excitation energies of 20, 40, 60, 80, 
100, and 120 meV, respectively. k x and k y both vary from 
-7r to 7r. Note k x is of spin-antiparallel direction and k y is of 
spin-parallel direction. 



tures are consistent with the strong anisotropy of the 
spin wave excitations along the spin-parallel and spin- 
antiparallel directions due to the Fe vacancies. 

For TlFei.5Se2, we also find the quite similar features 
on the spin dynamical structure factor (see Fig. [SJ. 
Since the spin wave dispersions become further weaker 
along the spin-parallel direction in TlFei.5Se2 than in 
KFei.5Se2, as shown in Fig. 3, the diffraction peaks thus 
become line-shaped along k y at 60 meV rather than 100 
meV for TlFei.5Se2. Moreover, when uif increases from 
80 to 120 meV, the diffraction peaks become pear-shaped 
patterns. The large green areas near the right and left 
boundaries of the extended Brillouin zone at 100 meV 
come from the flatness of the highest optical branch. 

I 



V. CONCLUSION 

We have shown that we can analytically solve the spin 
wave excitations for a complex magnetic structure by us- 
ing the equation of motion method in the framework of 
the linearized spin wave theory, which is illustrated by 
studying the magnetism of the newly discovered interca- 
lated ternary iron-selenide AFei.5Se2 (A=K, Tl). We find 
that there are one acoustic branch (gapless Goldstone 
mode) and two gapful optical branches of the spin wave 
excitations with each in double degeneracy in AFei.5Se 2 . 
The phase boundary of AFei.sSe2 in the 4x2 collinear 
antiferromagnetic order is determined by examining the 
non-imaginary excitation frequency condition, which in- 
corporates the antiferromagnetic quantum fluctuations. 
We also derive the exchange couplings between the Fe 
moments based on the first-principles total energy cal- 
culations, so that we can calculate and discuss the spin 
dynamical structure factors in connection with neutron 
inelastic scattering experiment. By computing the an- 
tiferromagnetic quantum fluctuations, we find that the 
Fe spin is S = | in AFei.sSe2 and a very small spin- 
orientation anisotropy can remarkably suppress the an- 
tifcromagnetic quantum fluctuations. 



ACKNOWLEDGMENTS 

This work is supported by National Program for Basic 
Research of MOST of China (Grant No. 2011CBA00112) 
and National Natural Science Foundation of China 
(Grant Nos. 11190024 and 91121008). 



VI. APPENDIX I: COEFFICIENTS IN EQ. (5) 



16 Jf cos( 



2C 2 - 2D? - D 2 , 



c = 4J 4 b - 32 J-yf cos(6 x ) 2 + 64 J 2 4 cos(0 x ) 4 

-2C 2 D 2 + Df + 16 Ji fc J 2 cos(^)C k D 2 - U 2 b D x D 2 



AJlCl 



16 J 2 cos 



fCl + C 4 + U\ h D\ - 16J| cos( 
16J 2 2 cos(6y 2 DiD 2 - 2Clpl + 2D\D\, 



d = 2Jf b cos(46gc£ + 2 Jf b C£ + 16J$,J$ cos(0 x ) 2 cos(40„)C£ - 16 J 2 b J 2 2 cos{6 x ) 2 Cl + 32 J 2 4 cos{8 x ) 4 cos(4^)Cg 
-32 J\ b Jl cos{9 x ) 2 Cl + 32Jlcos(9 x ) 4 Cl - 16Jf b J 2 cos{8 x ) cos(46gC k Di + 16 Jf b J 2 cos{6 x )C li D l 
-64J u Jl cos(6y 3 cos(46» !/ )C k Di + 64J 16 Jf cos(6y 3 C k Di - 4J 1 4 b D 2 + 32Jf 6 Jf cos(6y 2 cos(46» y )D 2 



-64j|cos(6 
-4J 2 6 D 3 D 2 



) 4 D( + 16 Ji 6 J 2 cos(^)C 3 D 2 - 4 J 2 b C 2 DiD 2 



16Jf cos( 



16Jf cos( 



l DfD 2 - G 4 Dl + 2ClD\D\ - D 4 D 2 2l 



nYClDiD 2 - 16 Ji fc J 2 cos(fy C k D 2 D 2 



r 



where C k , D\, D 2 , 
the text. 



and 9 y are defined in Eq. ([2]) in 
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